FINAL REFORT 


¥ 


i 


% 


L' 1 . 1.1 TED TO: 


TNSTITl TION: 


National Aeronautics arid 
Space Administration 
Langley Research Center 
Hampton, Virginia 2306b 


Department of Physics and 
Engineering Studies 
Hampton Institute 
Hampton, Virginia 23668 


TITLE 01* GRANT: 


Local Effects of Partly- 
Cloudy Skies on Solar and 
Emitted Radiation 


GRANT NUMBER: 


PERIOD COVERED BY THIS REPORT: 


CO-PRINCIPAL 


INVESTIGATORS 


RESEARCH ASSOCIATES 


NAG 1-67 


August 1, I960 to 
July 31, 1983 


D. A. Whitney 
D. D. Venable 


T. J. Griffin 
J. R. Foreman 



(NASA-CR-173189) LOCAL EFFECTS OF 
PARTLY- CLOUDY SKIES ON SOLAR AND EMITTED 
RADIATIONS Final Report, 1 Aug. 1980 - 31 
Jul. 1983 (Hampton Inst.). 93 p 

HC A05/M.F A0 1 CSCL 04B G3/4 7 


N8 4— 16711 


Dnclas 

00541 


LOCAL EFFECTS OF PARTLY-CLCUDY SKIES ON 
SOLAR AND EMITTED RADIATIONS - FINAL REPORT 


D. A. Whitney 
T. J. Griffin 
J. R. Foreman 

Department of Physics and 
Engineering Studies 

ABSTRACT 

The third year of the proposed three year project has now elapsed. 

The computer automated data acquisition system for atmospheric emittance, 
and global solar, downwelled diffuse solar, and direct solar irradiances 
has been fully operational for about two and one-half years. Hourly- 
integrated global solar and atmospheric emitted radiances have been measured 
continuously from February 1981 to August 1983. Hour ly-inte grated diffuse 
solar and direct solar irradiances have been measured continuously from 
October 1981 to August 1983.- One-minute integrated data have been made 
available for each of these components from February 1982 to August 1983. 

Atmospheric aerosol and turbidity measurements for the period February 
1981 through July 1983 have been analyzed and the results are presented here. 

The correlation of global insolation with cloud cover fractions for the 
first complete year's data set was completed. A theoretical model was developed 
to parameterize the effects of local aerosols upon insolation received at the 
ground using satellite radiometric data and insolation measurements under clear 
sky conditions. A February data set, composed of one-minute integrated global 
insolation and direct solar irradiances, cloud cover fractions, meteorological 
data from nearby weather stations, and GOES East satellite radiometric data 
was collected to test the model and used to calculate the effects of local 


aerosols . 
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INTRODUCTION 


A solar energy measurement station was established at Hampton Institute, 
February 15, 1981. Routine hourly Integrated measurements were made of global, 
diffuse, and direct solar irradiances and of atmospheric emittance. After the 
data acquisition was computer automated in February 1982, one-minute Integrated 
radiometric data, as well as one-hour integrated data, were recorded. More 
detailed information about the measurement system is presented in Sections I 
and II of this report. Monthly averages for global, diffuse, and direct solar 
irradiance, atmospheric emittance, and atmospheric aerosol and turbidity 
parameters were calculated and are presented in Section III. 

Correlation of global insolation with cloud cover fractions were made 
using the ARL model and the results for the first complete year of data are 
presented in Section III . A parameterization method for estimation of the 
effect of aerosols upon insolation has been developed and a data set has been 
accumulated to test this method. Detailed information about the parameterization 
method is provided in Section IV. 


I. SOLAR RADIATION MEASUREMENT PROGRAM 


An observation platform for solar and atmospheric radiation measurements 
was established on the roof of Turner Hall (latitude 37.02% longitude 76.34% 
and elevation 24 meters) February 16, 1981. A radio tower and a smokestack 
are the only two obstructions greater than ten degrees above the horizon with 
all other obstructions less than five degrees above the horizon. Routine 
measurements were made of global solar irradiance, diffuse solar irradiance, 
direct solar irradiance, and atmospheric emittance. Information (Griffin, 1982) 
about the solar radiation rreasure^nt program was presented at the April 20-23, 1982 
meeting of the Virginia Academy of Science held in Blacksburg, Virginia and 
at the June 1-3, 1983 meeting of the American Solar Energy Society in Minneapolis, 
Minnesota (Whitney, 1983). The Abstract of the presentations by T. J. Griffin 
are attached as Appendices I and II . A surmary of the solar and atmospheric 
data available at the time of this report is provided in TABLE I. 


A. Instrumentati on 

Since detailed descriptions of the radiometric instrumentation are 
available in the first Annual Report (Whitney, 1981) for this grant, only a 
surmary of our measurement capabilities is presented here. The global solar 
irradiance on a horizontal surface was measured by an Eppley Precision 
Spectral fy rename ter (PSP) with a hemispherical WG7 Schott glass dome. Diffuse 
solar irradiance on a horizontal surface was measured by an Eppley PSP with a V,'G7 
dome that was mounted on an Eppley Solar Tracker and Occulting Disk System. 

Direct solar irradiaice was measured by an Eppley Normal Incidence Pyrheliometer 


(NIP) mounted on a Solar Tracker. An Eppley Hickey-Frieden Absolute Cavity 
Pyrheliometer was used to calibrate the NIP regularly. Atmospheric emittance 


was measured with an Eppley Precision Infrared Radiometer (PIR). In addition, 
turbidity (Volz, 1974) measurements were made approximately hourly for deal 
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sky conditions using a Volz Sunphotometer. 

The wavelength range of each of these Instruments is listed in Table XI. 
Information about the measurement frequency and about the time periods that 
insolation data are available is presented in TABLE I. 

B. Calibration 

Each radiometer was calibrated prior to delivery by the manufacturer. 
Subsequent calibrations have been performed by comparison to secondary or 
primary standards at Hampton Institute or at Eppley Laboratory. Calibration 
data are presented in TABLE III by listing each instrument, the date and site 
of each calibration and the calculated sensitivity factor. The NIP was 
calibrated by comparison to the Hickey-Frieden Absolute Cavity Pyrheliometer 
which is considered to be a primary standard. There were no adjustments in 
the NIP calibration factor since the calculated sensitivity showed no change 
within the sensor accuracy. The pyranometers were compared twice a year with 
one another by comparison of three— day integrated global irradiance totals af'ter 
side-by-side operation. Comparison between the pyranometers indicated a consistent 
difference of about one percent which was within the two percent accuracy of 
each instrument, thus no adjustment in either calibration factor was necessary. 

A list of equipment used in these calibrations is attached as Appendix III . 

The pyrgeometer (PIR) and the pyranometers were recalibrated at 'Eppley 
Laboratory twice during the three— year period of use. A change in calibration 
standards at Eppley Laboratory in October 1981 accounted for a 2.6 percent 
change in sensitivity factor for the PbP instruments, and required that data 
obtained previous to that date be adjusted to standardize our data set. 

Annual calibrations of both PSP's at Eppley Laboratory showed a sensor 
degradation in each instrument of approximately 0.1% per month. A linear 
correction factor between Eppley calibrations for each instrument has been 


¥ 
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calculated and used for the monthly average measurements presented In this 
report. Raw data stored on magnetic tape have not been corrected for the 
above changes in sensor sensitivity. Monthly correction factors used to 
adjust all of the raw data on tape are presented in TABLE IV for each 
radiance component. 

The recorder systems were calibrated every six months by using a stable 
millivolt source and by adjusting the integrator and strip chart recorder 
gains to obtain the proper readings. The electronic integiators exhibited 
extreme stability with the maximum adjustment required for any integrate 
being only 0.4%, while larger adjustments were occasionally required for the 
chart recorders. 

The meter on the Volz Sunphotometer v/as replaced on April 29, 1983 
immediately after it became inoperative. Telephone discussions with the 
manufacturer revealed that the change in meter would have a negligible effect 
on the calibration of the instrument. An attempt to verify this by checking 
for linearity of the sensor on a Langley Plot of meter readings versus air 
mass failed due to an inability to obtain enough clear sky measurements on a 
given day over a significant range in air mass. Atmospheric aerosol and 
turbidity data reported in this report have been calculated using the original 
calibration data. 

C. Meteorological Data 

Standard hourly meteorological observations taken at nearby Langley Air 
Force Base (LAFB) were picked up on a monthly basis from Detachment 7, Third 
Weather Squadron and are on file at Hampton Institute. These data included 
information about cloud height, cloud type, fractional sky cover, precipitation, 
sea level pressure, and temperature. These data were supplemented by whole— sky 
photographs and visual cloud observations at Hampton Institute. Additional 
meteorological, turbidity and ozone data were purchased from the National 


-3- 


Climatic Center archives in Asheville, North Carolina and the World Ozone Center, 
Environment Canada, Toronto, Canada for use in the data analysis. 

D. Cloud Cover Fractions 

Cloud cover composed of cumuliform clouds were selected for this study 
because of the nearly opaque optical properties, low cloud altitude, distinct 
boundaries, and frequency of occurence. Cloud cover fraction was defined for 
the local area as the ratio of cloudy area to the total area. There are three 
different sources of cloud cover fractions: 1) satellite data, 2) ground 

based photography, and 3) trained observers. Only a summary of these methods 
is provided here since details about the analysis methods are provided in the 
first Annual Report (Whitney, 1981). These cloud cover fractions were used in 
the ARL regression equation as discussed in Section III-C of this report. These 
methods are compared with one another in the first Annual Report. 

1. Satellite Derived Cloud Cover Fractions 

Black and white photoprints of visible imagery provided by GOES-EAST 
were selected on the bases of: 1) frequency (every one-half hour), 2) range 

of cloudiness (all fractions possible), and 3) convenience (our local 
geography was easily distinguishable in the prints). A distance scale was 
calculated from known landmarks on the photoprints for the east-west and 
north-south directions. Using this information an ellipse corresponding to a 
120 km radius horizon circle was drawn on a clear plastic overlay. The ellipse 
was further subdivided into grids that corresponded to 24 km x 24 km squares at 
the ground. The central grid was placed over the Hampton Institute measurement 
site and visual estimates of cloud cover fractions made for each grid. These 
fractions were then used to calculate the cloud cover fraction for the local 
area. Comparison of the fractions obtained by this method to the other two 
methods indicated that the photoprint method was the least reliable for cumuliform 
clouds and, thus, satellite photoprint analysis was only used during the first 
year of this program. 


2, Ground Based Photograph-Derived Cloud Cover Fractions 

A whole-sky photograhlc system was constructed using a 35 mm SLR camera 
body, an Aetna fish-eye adaptor lens, and a camera mount directed toward 
zenith. The camera system was calibrated by aiming at a large flat surface 
(a classroom wall) and carefully measuring the radial distance on the 
photograph for each known angular position. A linear relationship was 
observed between angular position and radial distance in the photograph as 
indicated in Figure 1 from the center out to 85 degrees. 

The 85.0 degree field of view about the zenith defines the local area 
as a circle of radius from 2 to 50 km depending on cloud altitude. Over 
this limited field of view the atmosphere can be treated as being flat and 
the cloud cover fraction can be calculated independent of cloud altitude. 

(See the first Annual Report for details). An analysis grid for ground 
based photographs and slides was developed using concentric* rings and radial 
sectors. During the first year, black aid white photographs were enlarged 
to fit the analysis grid overlay and a cloud cover fraction determined for 
each grid by visual inspection. Afterwards color slides were directly 
projected onto the analysis grid in order to aid in distinction between dark 
cloud bottoms and clear sky, and to reduce both processing and analysis times. 

The whole-sky camera system was modified in order to allow computer 
activation of the camera by adding an electronic shutter and an autowind 
system. The camera system was mounted on the Eppley Shadow Band Stand in 
order to eliminate the need for frequent adjustment of the sunshade. A 
clock and date card were placed on the inside surface of the shadow band 
to document the time of day and date of each photograph. Whole-sky 
photographs (color slides) were taken every one-half hour on selected 
weekdays when clouds were present without precipitation, from August 1981 
through December 1982. 
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3. Visual Observations of Cloud Cover Fraction 
Visual estimates of cloud cover were made by trained observers at 
nearby Langley Air Force Base every hour. The observations included the 
cloud cover fraction in tenths and cloud type In code. These observations 
were screened to eliminate hours with predominantly transparent or semi- 
transparent clouds. Visual observations made at Harrpton Institute during 
the first year were used in order to help interpret the photoprint and 
photograph analysis results. 

E. Data Handling and Quality Control 

A data storage procedure for the radiation data was devised to 
efficiently handle the data and ensure quality control . The Integrated 
radiometric data and times were initially stored on a Tektronix 4051 
microcomputer's internal magnetic tape unit. Then on a monthly basis 
these raw data files were transferred via computer hookup from the Tektronix 
4051 to the PDP 11/34 minicomputer where the data were permanently stored 
on 1600 bpi magnetic tape. All preliminary processing was done on the 
Tektronix 4051, while data analysis routines and application programs were 
performed on the PDP 11/34 system. The data were examined for errors by a 
computer program that located the gaps in the data and identified places 
to be investigated and corrected. 

The automated data acquisition system for global and direct solar 
irradiances was fully operational from February 1, 1982 through July 31, 1983. 
During that time integrated radiometric data for one-minute intervals were 
obtained directly from the Eppley Integrators and stored on magnetic tape. 
Specific data handling procedures and quality control for these data are 
discussed in more detail in the next section of this report. 
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Prior to the installation of the automated data acquisition system, only 
one-hour integrated measurements were recorded on magnetic tape. Printed data 
from the integrators were scanned on a dally basis for missing or problem data 
and incorrect timing caused by power failures or other electrical and mechanical 
malfunctions. Missing data were supplemented by the strip chart record when 
available. Approximately once a week these data were manually entered into the 
Tektronix 4051 microcomputer, Inspected for operator errors, and then transferred 
to the PDP 11/34 data storage tape. 


I 
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IX, AUTOMATED DATA ACQUISITION SYSTEM 


One of the major objectives of the second year of this research project 
was to develop and implement tin automated data acquisition system. This 
system reduced the manhour requirements, and the number of operator induced 
errors, involved v,ith data transferral from the integrator-recording system 
to final permanent magnetic tape storage. The hardware requirements, 
operating procedures, features, and performance of the system were provided 
in the second Annual Report i Whitney, 1982) ond are surmarized below. 

A. System Overview 


Ihe automated data acquisition system connected the radiometric sensor- 


integrator recording system with the microcomputer I/O capabilities via on 
interface box. The microcomputer read the integrated radiometric values 
from the BCD interface panels at a preset time Interval and recorded these 


data on the internal magnetic tape unit. Whole-sky photographs were triggered 
by computer command at preset times and a record of time and picture number 
were recorded on the magnetic tape. Integrated atmospheric omittance, and 


global, direct, 


and diffuse 


solar irradianees were recorded at one-minute 


intervals from 0400 EOT to 2000 EOT and at ten-minute intervals for the rest 


of the night. The data for each month were transferred to separate 1600 bpi 
magnetic tapes for permanent storage. 

B. System Description 


The data acquisition system was composed of three separate subsystems? 

(1) the radiometric sensor-integrator system; 

(2) the integrated data sampling and recording system; and 

(3) the data storage system. 


A flow chart of the data acquisition system is presented in Figure 2. 
The radiometric sensor-integrator system is located at the top half of the 
chart and includes the following five components: 


-a- 


I# 


I * 


U) the Eppley radiometers; 

(2) the Epplcy Integrators that sunmed the instantaneous readings from 
the radiometers; 

(3) the digltec printers that provided a hard copy printout of one-hour 
radiance values; 


(4) the X-y strip chart recorders that provided a hard copy of Instance- 
ous radiances; and, 

(5) the camera system that took the whole-sky photographs. 

The integrated data sampling and recording system located at the lower 
left-hand side of the chart includes: 


(1) the integrator signals that were provided at the BCD interface on 
the Eppley Integrators; 


(2) <-he interface box that centralized the data for computer access; 
(S' WB HOMs that interfaced the BCD data to the Tektronix 4051; 
t4) the Tektronix 4051 microcomputer that read the ROMs and activated 
the camera photocounter ; and, 


is; the magnetic tape cartridge where the computer stored the data. 

The data storage system is located at the lower right-hand side of the 
chart and consists of three components: 


(1) another Tektronix 4051 microcomputer with a R.S232 interface that 
acted as a link between the PDF minicomputer and the internal tape 
of the Tektronix computer; 

(2) the PDF 11/34 minicomputer system that read the data from the 
Tektronix 4051; and 


(3) the 1600 bpl, 
data. 


track magnetic tape for permanent storage of the 
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The hardware used In the Integrated data acquisition and storage system 
are listed in Appendix XV along with two other devices that were used in the 
development and testing of the interface box and software programs. Detailed 
Information about the operation of these devices and about the construction 
of the interface box is kept in a docunentation file in the Solar Energy 
Measurement Laboratory . 

0. Operating Procedure 

The data acquisition system operated on the Tektronix 4051 microcomputer 
using the computer program SOLAUTO written In BASIC by D. D. Venable. The 
program worked by comparison of the time provided by the Real Time Clock ROM 
Pack with times that were calculated from fixed time intervals entered into 
the program for each sensor or for the camera system. A copy of the program 
listing is kept in the documentation file in the Solar Energy Measurement 
Laboratory. Basic operation of the data acquisition system was discussed in 
some detail in the second Annual Report . 

The radiometric data were recorded on the magnetic tape in 72 character 
records. The first two-digit flag (0//) of the data record indicated the type 
of data string using: 01 for radiometric data, 02 for a photograph record, 

03 for a user message, and 04 for a system message. The next characters 
indicated the day of the week (three letter abbreviation) separated by two 
spaces from the date in the font: day-month-year (DD-MMM-YY) . The date was 

represented by two digits each for the day and for the year, and by the three 
letter abbreviation for the month. The next set of numbers (HH:MM:SS) in the 
data string were the time (EOT) represented by two digits each for the hour, 
minute, and second, and the remaining characters wore either a message or 
data. A radiometric data record had the form: 

01, DDD DD-MMM-YY HH:MM:SB, CL// , NNNNN , DF// , NNNNN , 1 R// , NNNNN , DR// , NNNNN 


The radiometric data were recorded as integrated totals starting from zero 
at midnight using the two letter abbreviation for each radiation component, 
a single digit activation code (1 for on, 0 for off), and five digits for 
each reading (NNNNN) . Global solar irradiance was abbreviated by GL: 
diffuse solar by DF: Infrared (atmospheric emlttance) by IR: and direct 

solar by DR. 

D. System Performance 

The performance of the computer automated data acquisition system v/as 
measured by calculation of the amount of data lost in comparison with the 
data recovered on magnetic tape. The performance record of the previous 
data acquisition system v/as also considered since both systems had some of 
the same causes for loss of data. For example, diffuse solar radiometric 
data were not collected during calibration periods for either pyranometer in 
1982. Intercomparison of the pyranometers was made only for the horizontal 
global solar orientation. Electrical storms interfered with computer program 
execution and occasionally stopped data collection until the system was restored 
to normal operation. Severe storms reset the printer times and the integrator 
count values, an effect that destroyed the printer data until the system could 
be restored to normal operation. The strip chart recorders were used to 
retrieve most of the hourly integrated data lost by computer failure. Mechanical 
solar tracking failures sometimes caused the loss of the direct and diffuse 
solar data. 

The hourly integrated-data recovery record is presented in TABLE V for 
the full period of the insolation measurement program at Hampton Institute. 

TABLE VI contains one and ten minute integrated data recovery information 
for the first five months of automated data acquisition. The measurement 
interval for global and direct solar irradiances was set at one-minute 
starting at 0942 ESI February 1, 1982. One-minute integrated data sampling 
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for the diffuse solar and atmospheric emittance began at 2038 EST, 

March 1, 1982 when the second half ox ^he interface box was completed. On 
March 13, 1982 the sampling interval at night was changed to ten minutes 
(from 2000 EST to 0400 EST) in order to reduce computer tape storage 
requirements. One Tektronix Data Tape Cartridge is used to store approximately 
six and one-half { 6 ' 4 ) days of insolation data. 

Two parameters used to measure the performance of the automated data 
acquisition system were the average time period between failures and the 
average length of time lost for each error. These parameters were referred 
to as Meantime to Fail and Downtime per Error, respectively. The Mean Time 
to Fail was calculated by dividing the total possible number of data records 
by the number of failures and by the data record sampling rate (while the 
sampling rate was constant). The Downtime per Error was calculated by 
dividing the number of data records lost by the number of errors and by the 
sampling rate. The results of these measures of performance were reported 
(Blakey, 1982) at the 39th Joint Annual Meeting of the National Institute of 
Science and Beta Kappa Chi Scientific Society held in Washington, DC, 

March 17-20, 1982. Rody Blakey, an undergraduate assistant on this project, 
used the data obtained during the first month of computer automated data 
acquisition to calculate the Meantime to Fail and Downtime per Error. The 
results of his analysis were an average of 8.7 days between failures and an 
average of 6.7 hours downtime per error. The large downtime per error was 
caused by the two nighttime failures that stopped program execution for 
several hours until the system could be returned to normal the following 
working day. (Mr. Blakey received a Third Place Award for his paper in the 
Mathematics and Computer Science Section of the Meeting and one of two general 

awards given.) A copy of the abstract of his paper f "tached to this report 
as Appendix V . 
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HI . RADIOMETRIC AND METEOROLQG I CAL DATA ANALYSIS 


Several types of data analysis have been completed using tine radiometric 
and meteorological data collected at Hampton institute. The first type of 
analysis completed was the calculation of average values for the various 
measurements. A second type of analysis v.as required to treat the raw 
turbidity measurements in order to obtain useful parameters such as optical 
depths and precipi table water. Correlation of data with empirical formulae 
v.as a third type of analysis performed. Analyses involving satellite-derived 
radiometric data are discussed in Section IV as part, of tine paramo tori nation 
method. The results of tine other antilysis methods used are discussed below. 

A. Radiometric Data Presentation 

Diurnal variability in the various solar insolation components as observed 
by comparison of the average hourly values for afternoons with values for 


mornings. Average daily global, diffuse, and direct solar irradiances aid 
atmospheric omittance are listed ;r. .Aril. \ 1- . Average hourly global, diffuse, 
and direct solar irradiances and atmospheric omittance are plotted for each 
month from July 19el through June 1 9J3 in Hgure 3 tlirougli Figure 9. Diurnal 
variability is indicated by the lack of symmetry in these graphs. Seasonal 
variability can be seen by plotting global solar irradianee for clear sky 
days selected from each season such as is dene in Figure io. Flots can also 
be made for values averaged over shorter time intervals vone minute to sixty 
minutes' 1 . 
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B. Atmospheric Aerosol Extinction and Turbidity Data Results 

Individual measurements obtained with the Volz Sunphotometer (serial number 

492) were used to calculate aerosol optical depths and turbidity parameters 

using the formulation outlined in the second Annual Report and by Volz (1974). 

Monthly average aerosol extinction coefficients at 380 nm, 500 nm and 875 nm, 

along with turbidity coefficient 3 o and wavelength exponent <* o are listed in 

TABLE VIII for the measurement period-March 1981 through June 1983. The number 

of days on which measurements were obtained for each month are listed also. 

o .. 

The annual pattern of the Angstrom turbidity coefficient, $ 0 is graphically 
illustrated in Figure 11. As shown in the figure, 3 Q values peaked in the 
midsunmer. A high value of represents a high aerosol concentration in the 
local air mass. Atmospheric aerosols include dust, smoke, sea salts and other 
suspended particles. The aerosol concentration decreased in the autumn through 
the sprang season as represented by the 1981 through May 1982 data on Figure 11, 
however, in November 1982, values sharply increased instead of decreasing 
as expected. This increase is due to the effect of volcanic particles released 
into the stratosphere from the eruptions of El Chichon in April 1982. The 
global spreading of this volcanic dust cloud appears to have reached the 
Hampton Institute study region (37 °N, 76 °W) during the month of November 1982. 
A slight reduction in aerosol concentration occurred in December 1982 and in 
January 1983, but the normal annual overall reduction in aerosol concentration 
indicated by the September 1981 to May 1982 data did not materialize. In 1983 
the monthly average 3 q values generally increased through June. Aerosol 
extinction data at 380 nm, 500 nm and 875 nm showed a pattern similar to that 
illustrated by 3 Q throughout the measurement period. 

The ^ngstrbm exponent, a o , is a measure of aerosol size distribution. 

A large a Q means that small particles dominate the aerosol population, while 
a small a Q means that large particles dominate, ^ngstrbm ^states that a typical 
t in 1961 
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value of « 0 Is 1.3. At Hampton Institute the daily average value of „ o ranges 
from 1.73 to - 1.39 and the monthly average data do not follow any consistent 

annual pattern. The equation used to calculate <x is* 

o 

a o = 1.3 - 4.07 log (6-875 / 13-500), 

where 0-875 Is the turbidity coefficient at 875 nm, and S-SCO is the turbidity 
coefficient at 5CO nm. From this equation it is apparent that « o is determined 
by this beta ratio and that » 0 is equal to 1.3 only when the 0-875 / 0-soo 
ratio is unity. Negative values of a o are obtained when the beta ratio is 
greater than 2.087 which occurs when the local aerosol concentration measured 
at 875 nm is much greater than that measured "at 500 nm. 

C ' and Diffuse/Direct Relationsh i p with Atmospheric Turbidity 

Monthly mean solar irradiance values listed in TABLE VII were used to 
calculate the average daily Global to Direct ratio for each month during 1982. 

The annual variation followed a pattern similar to that of the turbidity 
coefficient In Figure 12 both the Global to Direct ratio and turbidity 

coefficient 0 Q are plotted versus time of year using the 1982 data. When the 
local aerosol concentration increased, the direct sola,- irradiance contribution 
to global Irradiance decreased and the Global to Direct ratio increased. 

Under clear sky conditions, direct solar radiation In the earth's 
atmosphere is only affected by scattering and absorption due to aerosol 
particles and atmospheric gases such as water vapor and ozone. If the concentration 
of aerosols increases, the amount of radiation scattered out of the direct beam 
to form diffuse radiation also increases. The relationship between the Diffuse 
to Direct ratio and aerosol concentrations, indicated by the coefficient „ 
was investigated by plotting the average hourly diffuse to direct irradiance 
ratio and versus time of day. On certain days, as illustrated for May 18, 

1983 in Figure 13, the two parameters were nearly equal and followed a very 
similar variation throughout the day. 
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For the year 1982, 166 individual turbidity coefficient a measurements 

o 

were plotted versus the diffuse to direct irradiance ratio calculated for the 

hour of each turbidity measurement. As shown in Figure 14 values of 3 varied 

o 

by as much as a factor of two for a given Diffuse to Direct ratio value at the 
low end of the graph. A least squares linear regression was completed to 
obtain the relationship: 

B o « 0.6302 ( DF / DR ) + 0.074, where, 

DF/DR Is the ratio of the average hourly solar diffuse to direct solar irradiance 
values. The line represented by the equation above is also plotted on Figure 14. 
The correlation factor of 0.919 obtained indicated a close linear dependence 
between B o and the average hourly diffuse to direct irradiance ratio inspite 
of the large variation discussed above. 

y oefficient 3 Q data were plotted versus average ten-minute 
diffuse to direct irradiance date, for each month. The linear regression analysis 
for each month resulted in a wide range of agreement with the assumed linear 
relationship. The November. 1982 data set of five measurements was used to 
calculate a linear correlation of 0.989. However, the May 1983 data set of 
44 measurements produced a correlation coefficient of only 0.062, implying 
that for this month a linear relationship between 3 Q and the diffuse to direct 
irradiance ratio does not exist. These results suggest that under the proper 
conditions the turbidity coefficient 3 q is linearly related to the ratio of 

the diffuse solar irradiance to the direct solar irradiance, but do not identify 
the other important variables. 

D. Correlations with Cloud Cover Fraction 

An empirical model (NOAA, 1979) developed by the Air Resources Laboratory 
(ARL) of the National Oceanic and Atmospheric Administration was selected for 
correlation of the global solar irradiance data with cloud cover fractions. 

The two equations that relate global insolation to solar zenith angle and 
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opaque uloud cover fraction are: 

SRC - A q + cos ZA + A 2 cos 2 ZA + A 3 cos 3 ZA (1) 

and; 

SR - SRC i^Bq + B 1 OPQ + B 2 OPQ 2 + Bg OPQ 3 ^ 

+ b 4 RN) . 

SRC Is the solar radiation hourly value for clear sky conditions. SR is the 
solar radiation hourly value for cloudy sky conditions. ZA is the zenith angle 
at the midpoint of each one-hour interval. OPQ is the average opaque cloud 
cover fraction, RN is a rain term that is equal to one if some form of 
precipitation is reported, otherwise it is zero. 

The coefficients for clear sky conditions were calculated separately for 
mornings and for afternoons each month of the year in order to partially 
account for diurnal and seasonal variations in atmospheric turbidity, water 
vapor, and other such factors. The first and last partial hours of daylight 
were not included in the regression calculation. The coefficients for the 
second equation were calculated for mornings and afternoons combined using the 
data for the first full year of insolation measurements. 

1. Clear Sky ARL Model Results 

The clear sky data had to be analyzed before the cloudy sky data could 
be normalized by the expected clear sky values. Determination of clear sky 
hours was made using the Langley Air Force Base (LAFB) cloud cover observation 
data set. The number of totally clear sky hours was insufficient for meaningful 
determination of the coefficients of equation 1 for most months and, therefore, 
these data were supplemented with "nearly clear sky" data for which the strip 
chart trace showed no indication of clouds and for which the cloud cover 
fraction was less than two-tenths. These nearly clear sky data were selected 
from hours which were coded as clear at either the beginning or the end of the 
hour, or were coded as having low fractions of transparent, or semitransparent 
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clouds. These data were added in order to Increase the number of points in 

the data set, and also to extend the meaningful range of the curve fit to 

zenith angles for which totally clear sky data were not available. This was 

conoistent with the original use of the ARL model in the rehabilitation of 
SOLMET data. 

The results of the application of this regression formula to our clear 
sky data are presented in TABLE IX for each month and in TABLE X for the first 
full year of global insolation measurements. The clear sky data are plotted 
in Figure 15 to show the extent of the agreement of the data with equation 1 
for the one-year data set from March 1, 1981 through February 28, 1982, Most 
of the scatter in the data was caused by seasonal variations of the various 
atmospheric constituents. The individual regression coefficients for each month 
have relatively large probable statistical errors associated with them and cannot 
be compared easily with other coefficients for a different month or measurement 
site. The coefficient Aj, is the most accurate term and, for the one-year data 
set, .it has only about a two percent probable error while even the algebraic 
sign of A 2 and A 3 is in doubt. The relative accuracy of the fit is demonstrated 
in the figure and, also, by comparison of the standard deviation to the data 

which gives a five to ten percent uncertainty in the midday insolation values. 

2. Cloudy Sky ARL Model Results 

The regression coefficients for clear sky mornings and for clear sky 
afternoons were used in equation 1 to determine the expected clear sky 
irradiance for each one-hour interval. The cloudy sky data were then normalized 
by the expected clear sky values and fit to equation 2 by using a nonlinear 
least-squares method. The cloud cover fractions were obtained from three 
different sources: 1) visual observations by trained observers (provided by 

the 3D Weather Squadron at nearby Langley Air Force Base-LAFB) ; 2) analysis 

of ground-based whole-sky photographs; and 3) analysis of GOES-EAST satellite 
photoprints . 
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Cloud cover fractions were calculated from the LAFB observations by 
averaging the value obtained five minutes before the hour began with the value 
at five minutes before the end of each hour (Eastern Standard Time), Hours 
with predominantly cirrus cloud cover were not used and thus only 2,926 one-hour 
averages were selected for use during the first one-year period. The data and 
curve plotted in Figure 16 are for dry conditions: RN=0. In order to satisfy 

computer space requirements, this large number of data was further reduced by 
the calculation of the mean value; the mean value plus, and minus, one standard 
deviation for each 0.05 step in cloud cover fraction from clear sky-0 to overcast 
sky-1 (21 steps plus overcast with precipitation times three values each for a 
total of 66 values). The results of the analysis of these LAFB derived fractions, 
which are presented in TABLE X and in Figure 16, reflect the use of these 66 
values along with 26 values of partly cloudy skies with precipitation during the 
one-hour period. 

The cloud cover fractions obtained from whole-sky photographs were plotted 
in Figure 17 for the one-year period. Most of the fractions were averages of 
the results of the analysis of two or three photographs and differ from the 
LAFB fractions in several respects. Only the first 70-75 degree field of view 
about the zenith in each photograph was used to measure cloud cover fraction 
as opposed to the standard 90 degree visual observation. Thin cloud cover and 
high clouds were not weighed heavily in the photographic analysis and thus this 
method provided a better measure of opaque cloud cover fraction. The curve 
obtained from the LAFB data is drawn in Figure 17 for comparison with the 
photograph derived fractions. 

A number of cloud cover fractions were obtained by analysis of GOES-EAST, 
black and white, visual image, photoprints of the local area by using an overlay 
grid. The results are plotted in Figure 18 along with the LAFB curve. 

Photoprint derived fractions less than 0.25 appear to fit better than data 
from the other methods. Very small, low lying, clouds can dominate a ground- 
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based photo or visual (JlK.vtvat.ion, but may not, it vlribp in the satellite 
photoprint because of tin 1 limited resolutivn. Very .mall thin clouds ?re 
unlikely to shade the pyrca lometer long enough to seriously .affect the hourly 
insolation value. Thicker clouds with small tmrmi t taneu and large reflectance 
can appeal'' larger In the photoprint. This thickness effect shifts cloud cover 
fractions derived from photoprint analysis relative u» ground-based fractions and 
improves the fit. All three plots shew (,:• nsidci'at le deviations from the curve 
for individual points, but none of the points obtained from photographs or 
photoprints are significantly above the clear sky value in Figure 17 or 18, 
while several points in Figure 10 arc? over twenty percent above the expected 
clear sky value, 

The partly cloudy hours with precipitation data provide an interesting 
test of the treatment of precipitation in equation These data tire plotted 
in Figure 19 with the overcast data represented by a mean value plus and minus 
one standard deviation. The lower curve represents the same coefficients as 
the other plots, but with KIf-1 . The upper curve uses the same coefficients but 
changes the precipitation term in equation to u* 4 »01-vHN) which quarantees 
no effect at OPQ = 0 (clear sky) ;.md reduces tu equation r for the normal 
overcast (OPQ = 1) precipitation condition. 
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The U»;on(tlofa-tji|»iWo<a 1 *v,l„ul Ly j. *. iomton, 
for use with our data and this method in below in detail. The cion 

sky part of this model wan applied to a one-month, meteorological and radio- 
metric data net m part of Mp. Rwn-n proposed u-etoral dissertation in 
Atmospheric iWu.v at the U,iv»T,;ity el' Kieidgan. herd data net was 
ased to test and to develop fine detailn in the model for entimatlon of 
aerosol absorption. The mum, of lehmary m ■ n.deeted an the tent 
monti, for application of Uiim moel.-l f or 1 1 „ • foiiowihg reasons: (l) t)K , 

Bround albedo wan expected to remain nearly constant throughout tto mot , th 
providing that tliere was no snow and that i.h, ,v was mu,. ln t|)l , 

amount 01 tee in the tidal banin; ta the aulomated data acquisition system 
provided short-,. I me interval lonc-minuhd integrated data; and (a) clear sky 

conditions wore sufficiently common to provide .. „ 

jjjiiV ut .1 data U«ii> tor eharac tort nation 

of ground albedo arid normal atmospheric ab.-orption and scatter Lrif.. 

A. Summary of the Model 

The method of estlnutiu, of global solar irraditnee on a horizontal 

surface is an extension of the short-waveh „g U , .rhergy hahneo equation 

developed by Ellis and Vender Hour { 1 t'/p) 

1 


‘hf! = (1 ~ a) (I 0 


r- "" 1 a) 


a 

I. 


hg iS thc horizont, ' J l >- Jobal wliortwivo irrndianec at, the earth's 
surface ; 

Is the local ground albedo; 

is the horizontal extraterrestriai solar irradimeo weighted by 
the spectral response of the satellite imaging device; 

IS the total Shortwave irradiaiico reflected to space t,y the 
earth's atmosphere arid , 


7 
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ft equals the tub'd ( nt and n*fp * 1 v >r 1 v.v iv' • radiation 

absorbed by the various components of Uu 1 at.mus{ here , 

The horizoi ital extraterrestrial insolation data were interpolated from 
measurements taken during the study pc Tied by *ui nbuulutt* cavity radiometer 
on board Nimbus 7 find reported by Hbi.* y, .* 1 . ai . two. The term l was 
calculated from the satellite measurd brightness »f the line eUm*nts and an 
appropriate bidirectional reflectance model for land and water surfaces 
developed by Rasehke, et. al , U‘>73). 

The portion of I, t due to f absorption by water vapor won estimated using 
total preclpi table water from nearby PAV/UlffjNbK mufcshrcmonts of the vertical 
dewpoint distribution. The absorption of critic and vat, op vapor was calculated 
from parameterization? ; developed by Laeis and Hanot 11 ( 1 974 ) . The absorption 
by permanent gases such as carbon dioxide find oxygon wore Cftlculated using a 
model presented by Burch, et. al . (lhfb) and by Yamamoto (loop). 

The absorption oi light by aerosols was computed as a residual from the 
measured four-minute global insolation data find the calculated absorption and 
reflection components » These calculated residuals were then parameterized to 
the column-mean relative humidity i, an found from a RAWINUQIJPK sounding. 

The weighting was done using a mean vertical distribution of aerosols in the 

low levels selected from iigurc one of shuttle anil ionn (lt)7t ) according to 

0 

the Angstrom turbidity parameter, beta. 

A fourth absorption portion of l ;i Is caused by clouds find was set to 
zero for this clear sky data set. 

B. Data Bets 

The meteorological data needed for this model were purchased from various 
sources. Most of the airways reports and weather data were purchased through 


ORIGINAL PAGE 13 
OF POOR QUALITY 


tho Environmental Data Cervices, national cj i untie c. nt*r, A::*! «»vi lie, Worth 
Carolina, Thu local study region in indicated in Hgurc 20 by the* dofctod 
linos find with tho Hampton institute me,'isureritent site indicated by an X. As 
noted in the legend, each meteorological site is represented by a separate 
number. The sources of data outside the study region were used to establish 
parameters at the boundaries of the study region .usd as- back-up data sots since 
several data sources did not report hourly* A summary of the data sets used 
in the analysis is given below in out line iVairi , 

1, Satellite Data 

a. Source: COEC-Fast visible and infrared digitised bright: less values . 

b. Resolution: 0.9 km x u.o urn visible, 3.7 km x 3,7 ton Infrared (at nadir). 

c. Data Array Coordinates: 


Patel 1 i to 
t Horn Inal ) 

Point line x element 
Center - gone x 749;* 

W Corner - 294 1 x 74 1 1 

PE Corner - , 'M4 l x 7 l b;-‘ 

HE Comer - 304b x 7 C i 2 

Cl V Conic r - 304e x /4‘ 1 


Coographical 

latitude x longitude 
37.01 o u fi x 7o . 338°W 
37 . f,87°H x 7b.8B4°W 
37 3 84°H x 7‘ .801°W 
3n.470°H x 7 l .79H°W 
3b.473°H x 7S.t c 8°W 


2 . 


d. Period: Clear fky Hours for February 7, h, 11, 13, It, 20 , 32, 23, 25 

and 28, 1982; Cloudy Cky Hours for* February 4, 13, 17, 18, 21 and 
24, 1982; 

Meteorological 
a. Insolation: 


i. Type: horizontal global - 3,3 m to 2 .H fJ m. 

11. Frequency: ten-minute totals centered on sate 11 Ite scan time (CMT). 

iii. Source: Kppley PCI' one-minut,e integrated data at Hampton institute, 

b. Humidity 


i. 


Hourly surface airways report 
altimeter set Lings from local 


; of b mpi-rat.uro, dew point, find 
airport;, Coast board, Air Force, 


Navy, 


and Army weather data sources 


for low level data. 


ORIGINAL PAGE [3 
OF POOR QUALITY 


ii. 


Radiosonde data frail Wallops to laid 


Dlerling Virginia, Capo 


I lab Leras aid Greensboro, North Carolina for upper level data. 

c. Cloud Parameters (not used in the clear sky analyses) 

1. Cloud fraction: Hampton Institute photographs aid satellite 

bri p^itneos readings . 

ii. Cloud type and height: Langley All’ Foret' Rase observations. 

Hi. Water content: radar reports (facsimile charts) from local 

airports . 

iv. Cloud top temperature: Infrared satellite data. 

d. Ozone and Aerosols 

i. Turbidity readings at Hampton Institute (or from a network of 

stations east of the Mississippi River, including, Rale] gh-Durham, 
North Carolina, for those eases, when turbidity reading, s were 


not tfikon at Hampton institute). 

ii. Dobson spectrophotometer readings of total ozone column over 
Nashville | Tennessee, Tal lahassee , Florida, Wallops Is laid, 
Virginia, and Washington, D.C. obtained from the World Ozone 
Data Center, Environment Canada, Toronto, Canada. 


Miscellaneous 

a. Ground albedo : calculated from clear sky satellite brightness 

values using an equation developed by Vender Haas and Fills U C )VL). 

b. Elevation: most points within the study region are treated as being 

at sea level out use of United dates Geological Gurvey Maps was 
made to ascertain elevations throughout the study region. 


c. 


Results 


A computer program "ABSOR.BAS", which solves the energy balance equation 
to estimate the absorption of wide-band short-wave radiation (hereafter called 
solar radiation) by aerosols was run for the Hampton Institute Solar Measure- 
ment Site (hereafter called SMS). The results for each case, Including all 
of the terms in the energy balance equation, are presented in Table XI. The 
fractional absorption (the absorption by aerosols, estimated from energy 
balance, and divided by the incident solar flux beneath the ozone layer), 
as well as the parameters used in each least-squares fit, are presented in 

Table XII for each case. These results are also presented in Figures 21 through 
23, where variables not explicitly named are assigned their mean values. 

One very salient feature of the absorptions by aerosols in both of the 
tables is that these are predominantly negative in value. This indicates a 
systematic overestimation of one or more of the terms of the energy balance 
equation; possibly the measured horizontal global insolation term (more on 
this later) or the gaseous absorption term, but most likely in the upwelllng 
reflection term owing to defects in the bi-directional reflectance model and 
in the image calibration. For this reason, the relative sun-satellite 
azimuth , Kasten's (1966) relative optical, air mass (hereafter called 
ROAM, which rs directly associated with the solar zenith angle) and the 
hours after 0000 GMT, January 1, 1982 were included in the least-squares- 
fit analyses. The sun-satellite azimuth and the ROAM were included to 
account for shortcomings in the bi-direotional reflectance model and the 
hours parameter was included to account for time-related "drifts" in the 
calibrations of the eight visual channels. 

Two of the tabulated cases were not used in any of the statistical 
analyses. Case No. 2 was omitted because of the presence of clouds 
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covering part of the 101 line x 101 element display sector allegedly 
containing the Hampton Institute SMS. Because of the high reflectivity of 
these clouds, the contrast of the non-cloudy portion of this sector was 
greatly reduced (the contrast of this display is made by dividing the range 
of brightness in the displayed sector into eight equal intervals and assigning 
a symbol to each interval) making the act of locating Hampton Institute very 
difficult and uncertain. Case No. 14 was omitted because, owing to the late 
hour of this image , the solar zenith angle was very la^ge requiring an 
unreasonable extrapolation of the bi-directional reflectance model. Also, 
because of the late hour, the visual image was very dark causing great 
difficulty in locating Hampton Institute. 

A third case, No. 19, however, could not be omitted on any such physical 
grounds, even though the fractional absorption for this case is lower than 
any of the others and much lower than any least-squares-fit formula prediction 
on this case, In fact, the inclusion or omission of this particular point 
made a great difference in the very nature and course of the least-squares- 
fit analysis and in the formula found from such an analysis. As a result, 
two sets of analyses, with and without case No. 19, were made. 

1. Case No. 19 included (two excluded cases): 

i. For turbidity related parameters excluded*: 

= 0.09322 64 - 0.000573405 0 

ss 

- 0.0635E6E f 


where i^ = fractional absorption by atmospheric aerosols; 

ii(-2) = two independent parameters with two cases excluded; 

®ss ~ relative sun-satellite azimuth (degrees) using the 

convention of Raschke et al. (1973) which defines 

*This was done to create a predictive formula which could be used in those 
large regions of the Earth's surface far from any Volz Sunphotometer 
observation sites. 
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this angle as the horizontal projection of the deflection of 
a photon from its incident direction (or the absolute value 
of 180 degrees minus the difference between the solar and 
viewing azimuths) rather than the absolute value of the 
simple difference between the solar and viewing azimuths; and, 
f = the mean relative humidity and is included as a "non- 
turbidity related parameter" since it can still be estimated 
without B, albeit less accurately, by weighing with a mean 
(B -independent) aerosol profile. 

The least squares fit analysis produced a multiple correlation 
coefficient, R = 0.44145 and an adjusted multiple correlation 
coefficient, R 1 = 0.3248G5. 

ii. For turbidity-related parameters included, the parameterization 
depends on the minimum acceptable value of p (-P[F - F ] (Os tie, 
1963; Bevington, 1969) chosen in finding the F-statistic: 
a. If p is chosen to be ^ 68%, only one parameter qualifies as 
an acceptable predictor and the result is: 

4" 2) = -0.0544005 + 0.138264 t 45 


where i(-2) = one independent parameter with two cases included 
and, 


t 4& = aerosol optical depth at x =0.45 pm. 


t 45 = m kp l ‘ 0 * 45 ) » 


where in = p'/po, 

m k = Raster's (1966) ROAM, 


p 1 = atmospheric surface pressure (kPa) at the point of 
interest (in this case, the Hampton Institute SMS), 
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and 

6 and a - Volz-$ngstrdJm turbidity coefficient and exponent, 
respectively. 

The analysis produced a multiple correlation coefficient, 

R * 0.541201 and an adjusted multiple correlation coefficient 
- 0.50565; 

b. If, however, p = 67%, almost all of the analysed predictors 

become valid in a stepwise least- squares process. In fact, 

the number of predictors used was limited to seven, not 

because of the computed F-statistic for predictors beyond 

the seventh, but because seven is the maximum number of 

independent variables which the multiple regression program 

("COSAP" statistical package) could handle: 

. vii(-2) 

X AA ~ "0.065407 + 0.367642 x - 0.000443137 e 

ss 

- 0 . 032541 7 m Rp 3 ( X^ («))“«+ o . 00863363 m kp 

- 0.157439 f + (9.87231 x 10 °) H - 0.00381737 a 
where 

\nax^ = Vi ’ ave l en gth of maximum aerosol attenuation of solar 
radiation found by setting the partial derivative with 
respect to X of the product of the Planck function and 3X“ a , 
equal to zero. After eliminating the zero and infinite 
roots, the resulting equation: 
exp (c'/X) = (o+ 5) X/( (a+ 5)X -c 1 ) 
where 

C = hc/'kT = 2.48067 ^um, 

h = Planck constant = 6.6256 x 10“ 27 erg s, 

c = speed of light in vacuo = 2.99725 x 10 10 cm/s, 
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k = Boltzmann constant = 1.36054 x 10" 10 erg/K, and 
T = mean effective solar radiative temperature = 5800 K 
(Glass tone, 1855), was numerically solved for A as a 
function of c, (very tricky, as there Is a very sharp — 
to +. singularity very close to the solution, especially 

at the larger values of a) , This function Is approximated 
to within 0.3% by: 

\nax^ = 0.0864642 - 0.00622149a + 1.75525/(a + 4,25) 


H Hours since 0000 GMT, January 1, 1982. 

The analysis produced a multiple correlation coefficient, 

R = 0.825926 and an adjusted multiple correlation coefficient 
R' = 0.71492. 


2. Case No. 19 excluded (total of three excluded cases): 


1 . 


For turbidity-related parameters excluded- 
, ii(-3) r 

M = 0. 3.32772 - 0.000797236 e _ (2.4183E x 10" 5 ) H 


M 

ss 

multiple correlation coefficient, R = 0.642801 and an adjusted 

multiple correlation coefficient, R< = 0,586281. 

ii. lor turbidity-related parameters included: 

. iii(-3) 

AA = 0.0485699 - 0.000488838 6 + 0.022923g a 

bb 

- 0.00646946 a 2 

The results were a multiple correlation coefficient, R = 0.84393 

and an adjusted multiple correlation coefficient, R. = 0.811331 

The results of all these least-s„ fits are presented in Ki^es 24 
through 28. 

rather surprising feature of all of these results is that wherever 
f appears, it has a negative coefficient whereas one would expect the 
opposite (e.g. MeWs, 197!; Hanel, 1972, 1976; Covert et al. , 1972; Pair 
andVorha, 1975; Fitzgerald, 1975; Fitzgerald et al. , i 982 ). 


One possible 


explanation of this is that there is a systematic overestimation of the 
absorption of solar radiation by water vapor, since the radlatively 
effective water vapor column on which the water vapor absorption estimation 
is based is strongly associated with Three possible sources of this are: 

1. The program "WATA1H.BAS" that estimates the radiatively effective 
water vapor column, w, from a sounding; 

2. The program "OZONAL.BAS" which fits a function of the form: 

z = c 1 + c 2 b + c 3 E + c 4 LE 

where z Is the quantity being fitted, L is the negative of the line 
number , E is the element number and c^,..., c^ are coefficients 
determined from the input data, using water vapor column observations 
from four upper-level stations just outside the study region. This 
function was used to obtain the radiatively effective water vapor 
column over Hampton Institute given Hampton Institute's line and 
element numbers in nominal GOES-east coordinates; and 

3. In the formula of Lacis and Hansen (1974) which estimated the 
absorption of solar radiation by water vapor given the local 
radiatively effective water vapor column and the relative optical 
air mass. 

One predictor which was applied to this analysis, but without success, 
was the formula modified from Hoyt (1978,1979) which purported to estimate 
the absorption of solar energy by aerosols: 
i^A - (1 ~ a5 s)d - g(3) mRp ) 

where 

w = albedo of single scattering by aerosols =0.95 and 
s 

g(e) = 0.937 - 1.044 (3 + 0.00575/(3 + 0.108) (which gives an 

acceptable fit to the tabulations of g(3) by Hoyt (1978,1979)). 

It is not known why this formula fared so poorly as a predictor of the 

absorption of solar radiation by aerosols. It may well be that the fault 
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lies with the modification of Hoyt's formula which was originally 1 = 

(1 - w ) fi(0) 1 • Numerous attempts since last October to contact Mr. Hoyt 

s 

for guidance on this matter have failed. Alternatively, this failure may be 
from errors in the estimation of absorption by aerosols from energy balance. 


V. SUMMARY AND COICU GIONb 
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Energy balance is a superb, fundamentally sound method, both for 
estimating the insolation at the Earth's surface and for investigating 
residual effects, such as the absorption of light by clouds or aerosols find 
errors resulting from defects in the bi-directional reflectance models aid 
in the calibration of the imaging device. In the case of this study, energy 
balance was used to investigate the absorption of solar radiation by aerosols 
and defects in the bi-directional reflectance model (Raschke et al . , 1973) 
and in the calibration of the eight visual channels of the GOES-east satellite 
(Norton et al . , 1980) . 

The results of this study are somewhat equivocal since they are 
excessively sensitive to the inclusion or exclusion of a single ease (No. 19). 
Obviously, a much larger data base is needed before some decisive conclusions 
may be reached, even for one location and during one month. 

Some facts about the results, however, are not so equivocal . For one 

thing, there is a decided negative dependence of the fractional absorption 

by aerosols on the relative sun-satellite relative azimuth angle , e , with 

ss 

or without case No. 19. The coefficient on e GC! is small in magnitude 
because e ss was given in degrees in this study. This strong dependence on 
0 ss i n< 3i- ca t es a systematic error in the bi-directional, reflectance model 
used in this study. This is also indicated by the dependence on the pressure- 
correlated relative optical air mass, m Rp (directly associated with tne solar 
zenith angle), in the predictive formula for p = 67% with case No. 19 and 
turbidity-related parameters Included. The dependence on the viewing zenith 
angle was not investigated in this study because of a limitation inherent in 
any view of a single Earth-surface point from a single geostationary satellite; 
namely, the viewing zenith angle varies very little, if at all. In fact, in 
this study, the range of this angle was less than two degrees. 


! 


Other strong dWnoieo found in thin study were on the computed 

optical depth at X = 0.4b /lm and on the VbtoAwutt. turt.Hlty exponent, 

especially on the square of «, in the predictive formulae including 
turbidity-related parameters, 

Other, weaker dependencies were found on: the aerosol prefile (selected 

by 3) weighted colimn mean relative htmidlty, f, a* aerosol optical depth 

at the “-dependent wavelength of maxim.™ absorption, m„ 6 ( x 

kp max' " ; 

and the number of hours, H, since 0)00 GMT, January 1, 198 2 which indicates 
a time dependence in the response (hence, in the calibration) of the satellite he 
visual channels. However, only a very weak dependence was found on a 
modification of the Hoyt (1978,1979) prediction, (1 - ®j (1 _ g(6) % } 

This parameter will, however, be retained as a regressive parameter to be 
investigated in future studies on a different or expanded data base. 

It has also been found that there be a large uncertainty, as much 

“ 3 ®’ ^ ^ 6Stlmatl ° n ° f the »«eeted spaceward by the Earth- 

atmosphere system. Of this, 20% may be due to problems in the calibration 

of the digitised visual Images from geostationary satellites (Muench, 1981; 

how this figurs of 20% was arrivpd o+- wi 1 1 k 

drived at will be discussed in Appendix IV. ) 

The remaining approximate 20 % (in a ovthamn.a,, > , 

^ Pythagorlan sum) is an educated guess 

on the uncertainty inherent in a bi-directional reflectance r^del and is 

comparable to the standard deviations *ich Davis and Cox ( 19 81) found in 

their own bi-directional reflectance model . 

another uncertainty arises from using four-minute totals from an 
integrator to obtain the measured mean horizontal global insolation at the 
Earth's surface. An integrator yields acceptable acctracy for totals of 
thirty minutes or more, but for shorter period totals, a different type 
of digitising device, or even an average of point measurements from the 
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analog trace, should be used, especially under low levels of daylight 
illumination (e,g., under cloud cover or near sunrise or sunset). This 
stands as a recommendation for future research as short-period averages are 
necessary in a high-resolution (in this case, the A-scale or 0,9km x 0.9 km, 
the highest nadir-point resolution available from an SMS/GOES satellite) 
study to ensure that cloud cover or illumination conditions do not 
appreciably change during the measurement. 

Another recorrmendation regards the bi-directional reflectance model. 

The model of Raschke et al. (1973) was chosen only because it was "tried 
and true", having been in the literature for a long time and used by 
numerous investigators in solar radiation and Earth radiation budget studies. 
The only other extensive model that was available at that time, Davis and 
Cox (1981), was rejected because the "bugeye" device used in that study 
sampled at only three nadir viewing angles (0°, 30° and 60°) whereas the 
nadir viewing angle in this study never strayed from the 40° - 43° range. 
Another model (Stowe et al., 1980) has been brought to our attention which 
may be superior to either previously identified model. 

Ideally, one should hire an aircraft and, using a photometric device 
as similar as possible in its spectral characteristics to that in tie SMS/ 
GOES VISSR, to make one's own survey of the angular (bi-directional) 
reflectance of the solar measurement site under as wide and as complete a 
range of solar zenith angles as possible during the study period. A 
further refinement would be to account for the angular reflectance of the 
overlying cloud-free atmosphere, including the changes of this angular 
reflectance with aerosol loading (e.g. , Braslau and Dave, 1973). Both of 
these refinements would greatly improve the accuracy and reliability of 
the estimates of the solar radiation reflected spaceward by the Earth- 
atmosphere system,, However, one considerable obstacle remains to complete 


accuracy In this estimation; the calibration uf ail oight bMf/UOb;; v.IDUR 
visual channels, Although numerous investigators ( 04 %, liauer and Lionesoh, 
1^75 > Hinton} Appendix 1 in Norton ot ul . , IQiiO; Mucnoh, 19U1) have been 
working on this very thorny problem, none have yet found a true solution. 

Another recommendation would be to use a solar measurement site as close 
as possible to a location where turbidity, upper level RAWiNSONPE and surface 
pressure (even Dobson spectrophotometer ozone) observations are made, 
thereby minimizing errors arising from extrapolating values of meteorological 
variables in space to the solar measurement site. Also, if the investigator 
has any control over such matters, the turbidity measurements should be made 
as close as possible in time to the expected satellite scan time (the time 
at which the VISSR on board the spinning satellite actually scans the 
measurement site as opposed to the nominal or actual image start time), 
thereby minimizing errors due to time extrapolation. Such extrapolations, 
both in time and space, can produce appreciable errors in the turbidity 
parameters, especially in an inhomogeneous atmosphere . 

A more complete study of the errors and their propagation in tire 
computations in this study will be made by Mr. Foreman in his Ph.D. thesis 
at the University of Michigan. Tills thesis is now in preparation. 
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TABLE I 


HAMPTON INSTITUTE SOLAR ENERGY MEASUREMENT SUMMARY 


Measurement 

Instrumentation 

Data Frequency 

Start 



SOLAR IRRADIANCES 





-EiJU. 

Global 

Eppley PSP with WG7 
clear glass dome 

One-hour integrated and 
continuous chart 

Feb. 17 

, 1981 

— 



Ore-minute Integrated* 

Feb. 1, 

1982 

JuL. 31, 1983 

Direct 

Eppley NIP with quartz 
glass and solar tracker 

One-hour integrated and 
continuous chart 

Oct. 1, 

1981 

— 



One-minute integrated* 

Feb. 1, 

1S82 

JuL. 31, 1983 

Diffuse 

Eppley PSP with WG7 
clear glass acme 

One-hour integrated and 
continuous chart 

Oct. 1, 

1981 

- 



Che-minute integrated* 

Mar. 1, 

1962 

July 31, 1S83 

ATMOSPHERIC EMITTANCE 

Eppley PIR 

Che— hour integrated and 

Mar. 18, 

1981 




continuous chart 





One-minute integrated* 

Mar. 1, 

1982 

JuL. 31, 1983 

ATMOSPHERIC PROPERTIES 





Aerosol Extinction 
@ 380 and 875 nm 

Volz Sunphotometer 

Approximately one hour 

Mar. 24, 

1961 


Turbidity @ 500 rai 


intervals for clear sky 




Precipltable Water 


Sd Sa ^ Ung rat ° - set *= on-Onute &», 0,00 EST 


TABLE l'I 

mdicmetric instrlmmtation WAVELENGTH CHARACTERISTICS 


Instrument 

Eppley Precision Spectral Pyranometer 

Eppley Precision Infrared Radiometer 

Eppley Normal Incidence Pyrhelicmeter 

Eppley Hickey-Frieden Absolute Cavitv 
Pyrheliometer y 

Additional Wavelength Ranges for the 
Pyranometers and Pyrheliometers 


Volz Sunpho tome ter 


Normal Observation Wavelength Range 
0,285 to 2,8 microns 


4.0 

to 

50,0 microns 

0,285 

to 

4.5 microns 

0.2 

to 

50 microns 

0.53 

to 

3.8 microns 

0.63 

to 

2.8 microns 

0.70 

to 

2.8 microns 

Center of Band 

- 

Halfwidth 

380 nm 

- 

11 nm 

500 nm 

- 

40 nm 

875 nm 

- 

17 nm 

940 nm 

— 

16 nm 


Radiometer 


TABLE III RADIOMETER CALIBRATION SUMMARY 


Calibration Information 


Date 


Conparison With: 


Performed By: 


Precis! * Spectral Pyranometer 


Sensitivity,, 
Factor ( uVm’l , r x ) 


Eppley PSP #20022F3 

10/1/80 

Standard References 


4/30/81 

PSP # 20613F3 


6/24/81 

PSP # 20613F3 


6/21/82 

PSP # 20613F3 


7/8/82 

Standard References 


9/6/82 

PSP # 20613F3 


7/7/83 

PSP # 22046F3 


7/18/83 

Standard References 

Eppley PSP # 20613F3 

2/28/81 

Standard References 


7/15/82 

PSP # 20022F3 


8/19/82 

Standard References 


7/7/83 

PSP # 22046F3 


9/9/83 

Standard References 
Standard References 

Eppley PSP # 22046F3 

10/13/82 

Standard References 

Normal Incidence Pyrheliometer 

Eppley NIP # 20254E6 

3/4/81 

Standard References 


5/12/82 

H-F Pyrheliometer 


5/5/83 

H-F Pyrheliometer 

Precision Infrared Radiometer 

Eppley PIR # 20078F3 

10/6/80 

Standard References 


7/7/82 

Standard References 


8/1/83 

Standard References 


Eppley Laboratory 

10.55 


Hampton Institute 

* 


Harp ton Institute 

* 


Hampton Institute 

* 


Eppley Laboratory 

lO.OOi 


Hanpton Institute 

* 


Hampton Institute 

* 


Eppley Laboratory 

9.82 


Eppley Laboratory 

11.10 

O O 
*n 2 

T3 £> 
O ^ 

Hanpton Institute 

10.62t 

Eppley Laboratory 

10.67 

Hanpton Institute 

■s 

O ;> 

Eppley Laboratory 


Z - 

‘-3 
, $ 

Eppley Laboratory 

10.52 

Eppley Laboratory 

10.47 


Eppley Laboratory 

9.21 


Hampton Institute 

** 


Hanpton Institute 



Eppley Laboratory 

4.95 


Eppley Laboratory 

4.86 


Eppley Laboratory 

4.90 



* '? iese sensitivity factors are unchanged within the ± 2% accuracy of the instruments 

The sensitivity factor was unchanged within the accuracy of the calibration instruments ( 1 5%) 
t Eppley Laboratory changed calibration standards October 1981 by 2.6%. ' ' 



HI 


* s 

TABLE IV RADIOMETRIC DATA CORRECTION FACTORS 


Month 

Year 

Global 

Direct 

Diffuse 

Atmospheric 



Irradience 

Irradiance 

Radiance 

Etalttance 

March 

1981 

1.033 

1.000 



April 

ti 

1.034 

(I 

_ 

1.0012 

May 

it 

1.035 

If 

- 

1.0024 

June 

ti 

1 .036 

II 

_ 

1.0036 

July 

ii 

1.029 

II 

- 

1.0049 

August 

it 

1.030 

M 

- 

1.0061 

September 

ii 

1.031 

II 

- 

1.0073 

October 

ii 

1.031 

II 

1.0413 

1.0085 

November 

ii 

1.032 

II 

1.0426 

1.0098 

December 

ii 

1.033 

M 

1.0439 

1.03 10 

January 

1982 

1.034 

II 

1.0452 

1.0123 

February 

n 

1.034 

II 

1.0464 

1.0135 

March 

ii 

1.035 

II 

1.0475 

1.0148 

April 

ii 

1.036 

II 

1_. 0488 

1.0160 

May 

n 

1.036 

II 

1.0499 

1.0173 

June 

ii 

1.037 

II 

1.0510 

1.0185 

July 

n 

1.038 

II 

1.0550 

1.0185 

August 

ii 

1.057 

II 

- 

1.0178 

September 

ii 

1.040 

If 

1,0582 

1.0170 

October 

n 

1.043 

If 

1.0598 

1.0162 

November 

ii 

1.045 

fl 

1.0614 

1.0155 

December 

ii 

1.0^8 

II 

1.0630 

1.0147 

January 

1983 

1.051 

II 

1.0646 

1.0140 

February 

ii 

1.013 

II 

1.0106 

0.9948 

March 

ii 

1.015 

II 

1.0121 

0.9940 

April 

ii 

1.018 

II 

1.0137 

0.9933 

May 

ii 

1.020 

II 

1.0152 

0.9926 

June 

ii 

1.023 

II 

1.0168 

0.9918 

July 

n 

1.025 

II 

0.9551* 

1.0183 


^Required for one-minute integrated data only for July 1983. 


TABLE V 


Year 


1981 


1982 


1983 


DATA RECOVERY RECORD; HOURLY INTEGRATED IRRADIANCES 


Month Number of Hour-Values Stored on Magnetic Tape 

Global Solar - Direct Solar - Diffuse Solar - Atmospheric Knitted 


February 

300 

— 


120 

March 

744 

— 

— 

664 

April 

720 

- 

_ 

718 

May 

744 

- 

— 

744 

June 

719 

- 

37 

720 

July 

741 

_ 

117 

741 

August 

740 

- 

0 

741 

September 

720 

13 

43 

720 

October 

743 

529 

712 

742 

November 

720 

615 

716 

720 

December 

742 

611 

732 

740 

January 

744 

676 

639 

722 

February 

672 

668 

672 

67? 

March 

743 

738 

744 

658 

April 

720 

692 

716 

720 

May 

740 

715 

721 

744 

June 

720 

720 

479 

512 

July 

744 

744 

279 

406 

August 

741 

741 

198 

742 

September 

720 

720 

567 

720 

October 

742 

742 

744 , 

744 

November 

720 

719 

720 

720 

December 

744 

740 

744 

744 

January 

742 

742 

738 

693 

February 

670 

671 

667 

665 

March 

744 

711 

742 

744 

April 

744 

699 

729 

744 

May 

744 

734 

744 

744 

June 

720 

690 

718 

711 

July 

744 

650 

744 

249 


TABLE VI 


DATA RECOVERY RECORD FOR AUTOMATED DATA ACQUISITION SYSTEM 
ONE AND TEN MINUTE INTEGRATED IRRADIANCES 


Month-Insol at ion 
Component 


Maximum Possible 
Number of 
Data Records 


Amount Recovered 
by Computer 
Number - Percent 


Number of Missing Data Records 
User System System 

Interrupt Crash Calibration 


Other 


February - 1982 


March 


April 


May 


June 


Global* 

Direct* 

Global 

Direct 

Diffuse** 

Infrared** 

Global 

Direct 

Diffuse 

Infrared 

Global 

Direct 

Diffuse 

Infrared 

Global 

Direct 

Diffuse 

Infrared 


40,320 

II 


36,648 

II 

It 

II 


30,240 

II 


31,248 

II 

II 

I! 


30,240 

II 

II 

II 


36,326 - 

90.1 

2,205 

1,190 


599 

35,163 - 

87.2 

3,358 

1,200 

- 

£99 

35,369 - 

96.5 

446 

678 


155 

35,381 - 

96.5 

433 

679 


ICC 

34,434 - 

94.0 

263 

679 


XvQ 

1 ' 7 ? 

26,470 - 

72,2 

257 

679 

- 

9,242 

29,673 ~ 

98.1 

61 

407 


99 

29,708 - 

98. 2 

26 

407 


99 

29,708 - 

98.2 

26 

407 


99 

29,708 - 

98.2 

26 

407 

- 

99 

29,519 - 

94.5 

17 

1,712 



29,519 - 

94.5 

17 

1,712 



29,519 - 

94.5 

17 

1,712 



29,519 - 

94.5 

17 

1,712 

- 

- 

28,232 - 

93.4 

181 

1,827 



28,232 - 

93.4 

181 

1,827 



16,112 - 

53.3 

181 

1,827 

12,120 


19,453 - 

64.3 

181 

1,827 

8,779 

- 

1, 1982 

to 2000 EST 

March 13, 

1982 and ten-minute readings 

at 


* One -minute readings only from 0942 Febru 
night after March 13. 

”S«terM™ 3 ? nly 2 ® 8 EST 1982 2000 ESr *** 13. 1982 and terminate readings at 
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TABLE VX1 


AVERAGE DAILY TOTAL IRRADIANCE SUMMARY 


Month 

Year 

Global 

Irradiance 

March 

1981 

(kJm"‘") 

15,610 

April 

ii 

19,814 

May 

ii 

19,919 

June 

ii 

22,604 

July 

n 

21,971 

August 

ii 

18,457 

September 

it 

18,378 

October 

n 

12,852 

November 

n 

9,292 

December 

n 

6,426 

AVERAGE - 

1981 

16,532 

January 

1982 

7,697 

February 

n 

9,965 

March 

n 

14,188 

April 

ii 

19,174 

May 

ii 

22,997 

June 

ii 

21,301 

July 

ii 

21,629 

August 

ii 

19,199 

September 

ii 

15,966 

October 

ii 

11,635 

November 

n 

8,831 

December 

ti 

6,131 

AVERAGE - 

1982 

14,893 

January 

1983 

8,017 

February 

ii 

10,539 

March 

ii 

13,150 

April 

ii 

17,605 

May 

ii 

22,624 

June 

ii 

24,421 

July 

ii 

25,027 


Diffuse 

Direct 

Irradiance 

Irradiance 

(kJnf 2 ) 

(kJnf 2 ) 

5,155 

- 

3,838 

— 

2 , 808 

— 

3,899 

10,552 

4,637 

10,634 

5,911 

13,219 

7,110 

15,095 

10,105 

17,708 

11,225 

12 , 848 

- 

14,087 

- 

13,450 

7,798 

12,866 

5,508 

12,013 

4,396 

10,300 

3,222 

7,538 

6,381 

12,526 

3,783 

9 , 990 

4,828 

10,958 

6,190 

10,163 

8,051 

12,254 

10,083 

16,330 

10,787 

17,395 

10,577 

19,681 


Atmospheric 

Emittance 

(kJirf^) 

19,242 

30,773 

33,420 

37,285 

37,746 

35,701 

33,264 

30,114 

27,288 

25,567 

30,940 

24,379 

26,600 

26,896 

28,458 

33,757 

36,169 

32 , 742 

31.086 

28.087 
26,532 
25,286 

29,090 

22,961 
22,945 
25,954 
26,902 
27,515 
31 , 648 
33,880 



TABLE VIII 


AVERAGE ATMOSPHERIC AEROSOL EXTINCTION AND TURBIDITY PARAMETERS 


Month/Year 

# of 

Aerosol Extinction Parameter 

- Tau 

Turbidity 

Parameters 


Days 

380 nm 

500 nm 

875 nm 

H 



Feb 1981 

3 

0.527 

0.127 

0.073 

0.064 

1.004 


Mar 

8 

0.353 

0.172 

0.856 

0.073 

1.176 


Apr 

6 

0.378 

0.172 

0.114 

0.105 

0.570 


May 

3 

0.455 

0.217 

0.144 

0.131 

0.643 


June 

5 

0.914 

0.547 

0.262 

0.220 

1.30 


July 

7 

0.871 

0.653 

0.353 

0.306 

0.942 


Aug 

7 

1.055 

0.667 

0.307 

0.256 

1.264 


Sept 

10 

0.612 

0.345 

0.173 

0.151 

0.888 


Oct 

3 

0.304 

0.107 

0.108 

0.109 

-0.004 


Nov 

2 

0.284 

0.105 

0.C88 

0.085 

0.291 


Dec 

2 

0.273 

0.129 

0.067 

0.058 

1.156 

o o 

Jan 1982 

5 

0.356 

0.146 

0.103 

0.096 

0.507 

-H IS 

Feb 

1 

0.327 

0.081 

0.105 

0.112 

-0.490 

T3 O 

Mar 

7 

0.477 

0.182 

0.146 

0.139 

0.379 

O 72. 
O ■£> 

Apr 

4 

0.332 

0.123 

0.133 

0.137 

-0. 252 

73 r* 

May 

2 

0.635 

0.254 

0.234 

0.230 

0.045 

lO -3 

June 

6 

0.940 

0.511 

0.323 

0.298 

0.596 

C 

63 

July 

5 

1.475 

0.925 

0.481 

0.414 

0.932 

r* pi 

Aug 

9 

1.402 

0.893 

0.482 

0.417 

0.943 

■**1 wrz*. 

<12 

Sept 

2 

0.654 

0.253 

0.158 

0.141 

0.829 

Oct 

4 

0.420 

0.171 

0.140 

0.135 

0.304 


Nov 

2 

0.625 

C.295 

0.198 

0.182 

0.523 


Dec 

2 

0.446 

0.235 

0.174 

0.162 

0.543 


Jan 1983 

5 

0.403 

0.194 

0.172 

0.167 

0.194 


Feb 

4 

0.556 

0.292 

0.214 

0.199 

G.5C8 


Mar 

2 

0.550 

0.233 

0.217 

0.214 

0.093 


Apr 

2 

0.670 

0.313 

0.270 

0.261 

0.244 


May 

7 

0.643 

0.151 

0.230 

0.258 

-0.868 


June 

4 

1.021 

0.436 

0.318 

0.298 

0.408 




TABLE IX 


ARL REGRESSION COEFFICIENTS FOR CLEAR SKY GLOBAL SOLAR IRRADIANCE 


Time of Day - Month 

Number 

of Hours 

Regression Coefficients CkJ 

*> 

m ~) 



Clear 

Only 

Nearly 

Clear 

A o 

4 

A 2 

A 3 


Mornings 

March, 1381 

10 

0 

- 212 

3360 

1256 

- 648 


April 

8 

11 

- 104 

2729 

2081 

- Pel 


May 

9 

*23 

43 

1422 

3992 

-1914 


June 

C 

8 

- 101 

2174 

1141 

36 C 


July 

6 

C 

- 94 

2423 

1022 

144 


August 

20 

C 

- 428 

4784 

-4068 

3337 


September 

35 

c 

- 86 

2268 

2797 

-1523 


October 

9 

0 

- 400 

5609 

-5652 

55C8 


November 

10 

10 

- 544 

2830 

1811 

- 482 


December 

1 

18 

- 17c 

3791 

18 

871 


January, 1982 

11 

1 

22 

1735 

6966 

-6289 


February 

1 

6 

- 144 

3424 

1102 

— 76 

O G 

Afternoons 

March, 1981 

14 

0 

- 187 

3856 

- 839 

1136 

-n X 

~0 c 

O 5. 

April 

5 

12 

- 50 

2419 

2740 

— jl cc 4 

O :r 

May 

0 

12 

173 

598 

6638 

-3834 

33 r* 

June 

5 

2 

- 7 

1184 

54C7 

-3143 

■O -.1 

July 

0 

6 

- 86 

1699 

3161 

-1206 

> 

August 

4 

2 

61 

£52 

5238 

-2351 

n 

September 

17 

4 

54 

756 

5738 

—2880 

63 

October 

1 

8 

- 32 

1735 

6185 

—5 404 


November 

8 

x 

- 270 

3949 

191 

148 


December 

c 

4 

- 486 

6163 

-6674 

7348 


January, 1982 

12 

2 

54 

936 

8766 

-6973 


February 

1 

8 

- 11 

2160 

3766 

-1897 



TABLE X 

ARL COEFFICIENTS FOR THE FIRST YEAR DATA SET 


(March 1, 

1981 through February 28, 

1982) 


Clear Sky Global 

No. of 

_2 

Regression Coefficients (kJ m ) 

Standard Deviation 

(Fit to Equation 1) 

Points 

A 0 

*1 

A 2 

A 3 

(in kJ m ) 

Mornings 







Clear Sky Only 

125 

-140 

3049 

1138 

- 616 

162 

Clear & Nearly Clear 

217 

-209 

3708 

- 302 

310 

158 

Afternoons 







Clear Sky Only 

72 

—252 

3920 

- 454 

554 

133 

Clear & Nearly Clear 

135 

-130 

2754 

2336 

-1328 

122 

Cloudy Sky Global 

No. of 


Regression Coefficients 

Standard Deviation 

(Fit to Equation 2) 

Points 

B 0 

B 1 

a 2 

B 3 B 4 


Reduced LAFB Data Set 

92 

1.004 

-0.353 

0.789 

-1.089 —0.215 

0.14 

Satellite Photoprints 

62 

1.017 

-0.468 

-0.182 

0.033 -0.253 

0.13 

Whole-Sky Photographs 

104 

0.955 

-0.615 

0.975 

— 0.85G -0.319 

0.14 
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TABLE XI 


Case 

Mo. 

Date 

(Feb. 

1982) 

Time 

(GMT) 

Surface 
Insolatii 
(Wh/irf } 

4 

7 

1407 

356.7 

C 

7 

1607 

620.4 

6 

7 

1807 

635.9 

7 

7 

2007 

418.8 

9 

8 

1406 

325.7 

lC- 

8 

1606 

604.9 

11 

8 

180C 

604.9 

13 

11 

2006 

418.8 

IE. 

13 

2207 

62.0 

17 

14 

It 06 

1 04 . 9 

19 

20 

1406 

403.3 

23 

23 

1406 

403.3 

24 

23 

1606 

666.9 

gsc 

2E 

1506 

CO J . 9 

26 

2E 

1637 

760.0 

27 

25 

1706 

77 c c 

23 

25 

1807 

744.5 

29 

25 

1907 

651.4 

30 

28 

1807 

729.0 

31 

28 

2007 

496.3 

32 

28 

2207 

108.6 


PARAMETERIZATION INPUT SUKMARV 


Extraterrestrial 

Opv.-elled 

Aerc-soi 

Gaseous 

Horizontal Inso- 
lation (Wh/nT) 

Reflectance 

(Wh/nT) 

Absorption 

(VJh/rcT ) 

Absorption 
(Wh/nr ) 

501.3 

91.19 

3.12 

81.24 

805. G 

155.19 

-22.47 

106.54 

833.4 

153.92 

-9.09 

107.80 

575.8 

108.95 

-1.71 

£6.17 

502.9 

113.73 

0.74 

90 . gG 

810.0 

' 52.18 

-16.04 

123.43 

840.2 

J 56.36 

4.42 

127.06 

603.3 

124. 84 

—7.83 

103.86 

140.8 

4c . 64 

— G . 7; 

140.77 

847.9 

157.12 

-4.34 

142.72 

569.4 

149.04 

—49.19 

IOC. 65 

582.4 

116. CO 

-7.94 

106.13 

897.5 

163.74 

-15 .45 

140.17 

783.0 

139.61 

—25 .!£. 

116.16 

946.9 

176.56 

-52.94 

129.25 

963.8 

202.34 

-74.60 

127.93 

943.6 

168.44 

-29.19 

124.48 

849.6 

145.20 

-5.36 

114.92 

960.3 

164.51 

-14.70 

144.79 

702.9 

133.16 

0.61 

115.94 

225.5 

62.04 

6.62 

55.66 
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TABLE XII PARAMETERI ZATI QTI OUTPUT AMD RESULTS 


Case 

No. 

Fractional 
Absorption 
by Aerosols 

Hasten ' s 
Relative 
Air Mass 

Optical 


. Relative 

Turbidity 

Average 

Depth x 
(0.45 /pm) 

Hoyt's 

Prediction 

Sun Satellite 
Azimuth (°) 

Alpha 

Beta 

Relative 

Humidity 

4 

0.00647 

2.850 

0.3760 

0.01067 

131.8 

2.309 

0.0595 

0.299 

5 

-0.02872 

1.781 

0.3616 

0.00676 

161.3 

2.303 

0.0575 

0.282 

6 

-0.01123 

1.719 

0.3177 

0.00622 

168.2 

2.210 

0.0544 

0.270 

7 

-0.00307 

2.477 

0.2029 

0.00882 

141.7 

1.631 

0.G552 

0.263 

g 

0.0G1E4 

2.830 

0.3326 

Q.G10G1 

131.4 

2.240 

0.0556 

0.331 

10 

-0.02294 

1.765 

0.3334 

G.GCC52 

160.8 

2.239 

0.055b 

0.318 

11 

0 . 00642 

1.699 

0.3334 

0.00631 

167.6 

2.236 

0.0559 

0.308 

13 

-0.01346 

2.354 

0.3553 

0.01486 

142.0 

1.491 

0.1060 

0.506 

16 

-0.00611 

9.266 

0.3696 

0.03686 

116.3 

1.590 

0. 1039 

G.3C2 

17 

-G.Q0£2fc 

1.675 

0.3819 

0.01041 

158.8 

1.675 

C. 101-3 

0.503 

19 

-0.08997 

2.457 

0.2326 

0.00457 

129.2 

2.708 

0.0268 

0.521 

23 

-0.01422 

2.414 

0.3063 

0.00561 

128.2 

2.750 

0.0341 

C.506 

24 

-0.01778 

1.571 

0.3095 

0.00356 

156.4 

2.823 

G..0325 

0.481 

25 

-0.03312 

1.814 

0.0832 

C. 01 141 

140.7 

-0.264 

0.102b 

0.513 

26 

-0.05746 

1.503 

0.0616 

0.01091 

164.2 

-0.817 

0.1183 

0.488 

27 

-C. 07944 

1.477 

0.0563 

0.01177 

172.4 

-1.062 

0.1314 

0.462 

28 

-0.03175 

1.509 

0.0894 

0.01075 

170.4 

-0.324 

C.1159 

0.436 

29 

-0.00649 

1.675 

0.0918 

G. 01100 

154.6 

-0.191 

0.1069 

0.410 

30 

-0.01570 

1.477 

0.3019 

0.00396 

170.3 

2.569 

0„03H8 

0.636 

31 

0.00090 

2.014 

0.3018 

0.00531 

139.2 

2.571 

0.0387 

0.315 

32 

0.04053 

6.083 

0.3019 

0.01438 

108.6 

2.573 

0.0387 

0.567 
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Whole Dky Photographic System Calibration Plot 

Automated Radiometric Data Acquisition System Plow Chart 

Average Hourly Global Irradianee find Atmospheric (110 Kmittanee 
for Each Month from March thru June 1981 


Average Hourly Global Irradianee and Atmospheric (IIO Emittanco 
for Each Month from July thru October 1981 

Average Hourly Global and Diffuse Col fir Irradionees find Atmospheric 
(IK) Kmittanee from Each Month from November 1981 thru February 1982 

Average Hourly Global, Diffuse md Direct Molar Xrradiances and 
Atmospheric (IR) Kmittanee for Each Month from March thru June 1982 

Average Hourly Global, Diffuse and Direct Molar Irradiances and 
Atmospheric (IR) Emittanco for Each Month from July thru October 1982 

Average Hourly Global, Diffuse and Direct Dolor Irrad lances and 
Atmospheric (IR) Kmittanee for Each Month from November 3982 thru 
February 1983 


Average Hourly Glob a, Diffuse find Direct Dolor i n-adiances and 
Atmospheric (IR) Kmittanee for Each Mon til from March thru June 1983 


Seasonal Variability of Clear-Dky Global Irradianee for Delected Days 

Monthly Average ftigstrdm Turbidity Coefficient B for the Period 
March 1981 thru June 1983 0 


Global to Direct Dolor Irradianee Ratio and Turbidity Coefficient 3 
for the 1982 Calendar Year 0 


Diffuse to Direct Dolor Irradianee Ratio and Turbidity Coefficient 
3 q for May 18, 1983 

Turbidity Coefficient 3 Versus Diffuse to Direct Molar Irradianee 
Ratio for the 1982 Calendar Year 

Clear Dky Global Insolation Versus Molar Fenith Anglo 

Cloudy to Clear Sky Global Insolation Ratio Versus Cloud Cover 
Fraction-Visual Estimates 

Cloudy to Clear Sky Global Insolation Ratio Versus Cloud Cover 
Fraction-Whole Dky Photograph Results 

Cloudy to Clear Sky Global Insolation Ratio Versus Cloud Cover 
Fraction-Satellite Photoprint Results 

Cloudy to Clear Sky Global Insolation Ratio Versus Cloud Cover 
Fraction-With Precipitation Present 

Hampton Institute Study Region and Local Meteorological Data Sources 


Fractional Aerosol Light Absorption m a Function of: (upper plot) 

Relative Sun-Satellite Azimuth, 9 ; find (lower plot) Raster's Relative 

Optical Air Mass, rn, . b,J 

Kp 
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< • Iraettonai Aem;ui J ,i a \l it An:.;< rrption or, a Inuct h jiy 4’5 (upp« r plot.) 
Aorocol-I rofilo Weighted Me; ip kiative Humidity, T'; ond, (lover plot) 
Voln-Angsti'Hm Turbidity < r if t *t , a, 

P3, fractional Aerosol Light Absorption an a Function of: (upper plot) 

Optical Depth at a Wavelength of 0,41 urn; arid, Clnv.fr plot) Hours, 
since (XXX; OMT, 1-H-HO 

,*4, bontter Diagran of Enietiunal bight. Absorpt ion by Aerosols Kstirnated 
from Energy balance Versus Predicted Absorption for Case 1,1, 

Ul. Scatter Diagram of fractional UgJIit Absorption by Aem:oUs Estimated 
from Energy balance Versus Predicted Absorption for C/rse 1,11 (a), 

c'6. Scatter Diagram of fractional bight Absorption by Aerosols Estimated 
from Energy balance Versus Predicted Absorption for Carre l.it (b), 

27, Scatter Diagram of fractional light Absorption by Aerosols Estimated 
from Energy balance Versus Predicted Absorption for Carre ?.i. 

28* Scatter Diagram of fractional Eight Absorpt ion by Aerosols Estimated 
from Energy balanct e-rsus Predicted Absorption for Case 2.11. 
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FIGURE 1. Whole-Sky Fhotographic System Calibration Plot. The distance 
of the image from the center of the photo is plotted versus the angular 
position. The error bars represent one standard deviation for each 
orientation for several points on each of tv:o calibration chctorrarhs . 
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FIGURE 5. Average Hourly Global and Diffuse Solar Irradiances 
and Atmospheric (IR) Emittance for Each Month from November 1981 
thru February 1982. Direct solar irradiance is plotted for 
January and February 1982. 
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FIGURE 6. Average Hourly Global, Diffuse and Direct Solar' 
Irradiances and Atmospheric (IR) Emittance for Each Month from 
March thru June 1982. 
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FIGURE 8. Average Hourly Global, Diffuse and Direct Solar 
Irradiances and Atmospheric (IR) Emittance for Each Month from 
November 1982 thru February 1983. 
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FIGURE 9. Average Hourly Global, Diffuse and Direct Solar 
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FIGURE 10. Seasonal Variability of Clear-Sky Global 
Irradiance for Selected Days. 
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FIGURE 12. Global to Direct Solar Irradiance Ratio and Turbidity 
Coefficient g for the 1982 Calendar Year. Irradiance values for 
the hour containing each turbidity measurement were used for each 
ratio . 
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SOLAR DIFFUSE TO DIRECT RADIATION RATIO & TURBIDITY B0 vs. TIME 
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1 HOUR AVERAGED DATA 

CIRCLE- ANGSTROM TURBIDITY PARAMETER, BETA 0 
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were used for each ratio. containing each turbidity measurement 
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SOLAR DIFFUSE/DIRECT RADIATION RATIO 

Trra^flnpp Coefficient 3 versus Diffuse to Direct Solar Irradiance Ratio. 

fS?^f ^urement are ccnpared to 8 
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FIGURE 15. Clear Sky Global Insolation versus Solar Zenith Angle. The data for 
turnings are plotted at the left and for afternoons at the rigit for the one year 
period March 1981 thru February 1982. Hie curve represents the result of the least 
squares fit to equation 1 for each data set. 
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«. ■ ■ OPAQUE CLOUD COVER FRACTION 

FIGURE 16. Cloudy to Clear Sky Global Insolation Ratio Versus Cloud 
Cover Fraction- Visual Estimates. Each hourly insolation value was 
normalized by the expected clear sky value,. Cloud cover fractions are 
averages of the hourly visual observations made at Langley Air Force Base 
The curve is the result of a least squares fit to equation 2. ' 
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‘ OPAQUE CLOUD • COVER • FRACTION 

FIGURE 17. Cloudy to Clear Sky Global Insolation Ratio Versus Cloud Cover Fraction 
Whole Sky Photographs. Each hourly insolation value was normalized by the expected 
clear sky value. The cloud cover fractions were obtained by analysis of black and 

white prints and of color slides of the local sky. The curve represents the result 
of the least squares fit to equation 2 using the LAFB data set. 
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OPAQUE CLOUD 'COVES j FRACTION 

•FIGURE 18. Cloudy to Clear Sky Global Insolation Ratio Versus Cloud 
Cover Fraction-Satellite Photoprints. Each hourly insolation value 

was normalized by the expected clear sky value. The cloud cover 
fractions were obtained by the analysis of GOES-East photoprints Th 
curve represents the result of the least squares fit to equation 2 
usinF the LAFB data set. 




CLOUDY/CLEAR SKY INSOLATION RATIO 


ORIGINAL PAGE IS 

OF POOR QUALITY 





U? 


ORIGINAL page fS 
OF POOR QUALITY 


FIGURE 20* Hampton Institute Study Region and Local Meteorological Data Sources 



Legend 

X- Hampton Institute 
l- Langley ' . F. B . 

*» Newport News, VA 

- Fort Eustis, VA 

- Milford Haven, VA 

5- Richmond , VA 

6- Sterling, VA 


7- Patuxent River ,MD 

8- Salisbury, MD 

9- Wallops Island, VA. 

10- Paramore Beach, VA 

11- Navel Air Norfolk, VA 

12- Norfolk, VA 

13- Oceana Navel Air, VA 


14- Cape Henry, VA 

15- Chesapeake Lighthouse 

16- Elizabeth City, NC 

17- Dare County, NC 

18- Oregon Inlet, NC 

19- Cape Hatteras, NC 

20- Rocky Mount , NC 
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LEGEND — PLOTTING SYMBOLS; 


» FOR TAU45 LESS THAN 0. 1 
+ FOR TAU45 BETWEEN 0, I AND 0 22 
X FOR TAU4S BETWEEN 0,20 AND 8 3 
0 FOR T All 45 GRFATFP THAN 0.3 


FIGURE 22. Fractional Aerosol Light Absorption as a Function of* funnpr> 
plot) o Aerosol-iWile Weighted Mean Relative Humidity, f* and (lower nlot) 
Volz~.Mngstrtim Turbidity Exponent, a. y ’ U r plot; 
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LEGEND — PLOTTING SYMBOLS •. 


FOR a. » 
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FIGURE 23. Fractional Aerosol Light Absorption as a Function of: (upper 

plot) Optical Depth at a Wavelength of 0.4b ym; and, (lower plot) Hours 
since 0000 GMT, January 8, 1982. 
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TWO-PARAMETER PREDICTED FRACTIONAL ABSORPTION 


FIGURE 2A. Scatter Diagram of Fractional Light Absorption by Aerosols 
Estimated from Energy Balance Versus Predicted Absorption for Case l.i. 
Tie fractional light absorption by aerosols was predicted by a least - 
squares-fit formula derived from the data set having two cases and all 
turbidity related parameters excluded. 
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ONE-PARAMETER PREDICTED FRACTIONAL ABSORPTION 


FIGURE 25. Scatter Diagram of Fractional Light Absorption by Aerosols 
Estimated. from Energy Balance Versus Predicted Absorption for Case l.ii(a) 
e fractional light absorption by aerosols was predicted by a least-squares* 
1 ? ?pj rnula derived from tde data set having two cases excluded and with 

P - > DO yu. 
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TWC-PARAMETER PREDICTED FRACTIONAL ABSORPTION 


FTGl'.RK 27 . Scatter Diagram of tract ional Light Absorption by Aerosols 
Estimated from Energy Balance Versus Predicted Absorption for Case H.i. 
The fractional l ight absorption by aerosols v/as predicted by a least- 
squares-fit formula derived from the data sot having three cases and 
all turbidity related parameters excluded. 
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THREE -PARAMETER PREDICTED FRACTIONAL ABSORPTION 


FIGURE 28. Scatter Diagram of Fractional Light Absorption by Aerosols 
Estimated from Energy Balance Versus Predicted Absorption for Case 2.ii. 
The fractional, light absorption by aerosols was predicted by a least- 
squares-fit i nriula derived from the data set having three cases excluded. 
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APPENDIX I 


.ABSTRACT - VIRGINIA ACADEMY OF SCIENCE PRESENTATION 


ENVIRONMENTAL SCIENCES 


DEVELOPMENT OF A SOLAR ENERGY MEASUREMENT LABORATORY FOR THE STUDY OF 
INSOLATION VARIATIONS AT HAMPTON , VIRGINIA. T. J. Griffin* . D. A. Whitney, 
and D. D. Venable. Dept, of Physics and Engineering Studies, Hampton Institute, 
Hampton, Virginia 23668. 

The purpose of this three-year study is to investigate the cloud 
dependence of Incident solar radiation at Hampton, Virginia. Solar 
irradlance at the Earth's surface is related to the extraterrestrial solar 
irradiance, to radiation absorbed and emitted by the atmosphere and clouds, 
and to radiation reflected by the Earth-atmosphere system. A ground-based 
measurement station has been established at Hampton Institute to monitor 
solar radiation, atmospheric emitted radiation, local cloud cover, and 
atmospheric turbidity. Continuous measurements of global, direct and diffuse 
solar radiation, and atmospheric infrared radiation are made and stored by 
computer. NOAA GOES-EAST satellite data are used to obtain albedo and cloud 
cover information. 

Interim analyses performed on the data include monthly averages of 
global insolation, infrared radiation, and atmospheric turbidity. Global 
insolation has been correlated with fractional cloud cover from March 1, 1981 
through February 1, 1982 using the ARL empirical model. 

(Supported by NASA grant No. NAG-1-87) 


APPENDIX II 


ABSTRACT - AMERICAN SOLAR ENERGY SOCIETY PRESENTATION 


SOLAR ENERGY MEASUREMENT PROGRAM AT HAMPTON, VIRGINIA 
D. A. Whitney, T. J. Griffin and D. D. Venable 


A global, diffuse and direct solar irradiance and atmospheric emlttance 
measurement program was initiated in February 1981 at Hampton Institute, 

Hampton, Virginia. Beginning March 1, 1982 the integrated irradiance and 
emlttance data were sampled on a one-minute basis and stored on magnetic tape 
by a microcomputer. Whole-sky photographs are used to document local cloud 
cover and are obtained on a regular basis. Atmospheric turbidity measurements 
are performed for clear-sky conditions with a Volz-type Sunpho tome ter. 

Several types of analysis have been performed with the radiometric data. 

Hourly global insolation has been correlated with opaque cloud cover fraction 
using the Air Resources Laboratory empirical model 1 . The cloud cover fractions 
were obtained from three different sources: 1) analysis of satellite photoprints; 

2) analysis of ground-based whole-sky photographs; and, 3) visual observations 
made by trained observers at nearby Langley Air Force Base2, Results of the 
comparisons for the first complete year of measurements will be presented. 

Mean hourly and daily total integrated irradiance will be presented for 
each month since February 1981. Atmospheric turbidity data have been analyzed 
in terms of the Angstrom turbidity parameters3 and aerosol optical depths at 
390 nm, 500 nm, and 875 nm. The results of the data analysis will be presented 
for the time period February 1981 through January 1983. 

This research was supported through the NASA Grant # NAG 1-87. 


1. NOAA, 1979: SOLMET Vol . 2., Final Report , Ashville, NC. USDOC/NOAA National 

Climatic Center. 

2. We gratefully acknowledge the assistance of the personnel of Detachment 7, 

3D Weather Squadron, Langley Air Force Base. 

3. Angstrom, A., 1961, "Techniques of Determining the Turbidity of the Atmosphere,' 
Tellus XIII , 2, pp. 214-222. 




APPENDIX III 

CALIBRATION INSTRUMENTATION USED AT HAMPTON INSTITUTE 


1. Calibration of Eppley NIP#20254E6 

Standard Sensor- ippley Model W-F Self-calibrating Cavity Pyrheliometer 
Serial Number 18752 

H-F Control Unit- Eppley Model 405, Serial Number 6621 

NIP Output Monitor- Keithley digital multimeter, model 179-20A, 

Serial Number 27764 


2. Calibration of Eppley Electronic Integrators 

Standard Millivolt Source- Honeywell Rubicon Potentiometer, Model #2730 

Serial Number "NASA-Langley 103291" 

Standard Source Monitor- Fluke Digital Voltmeter, Model #8300A 

Serial Number 307 
calibrated 9-10-81 

Eppley Integrator Voltage /frequency monitor- 

Hewlett-Packard Timer /Counter, Model 5327A 
Serial Number 1120A00231 
calibrated 12/83 

Integrator Amplifier Gain and analog output monitor- 

Keithley digital multimeter, Model 179-20A 
Serial Number 27764 



APPENDIX IV 


DATA ACQUISITION AND STORAGE HARDWARE 


ITEM DESCRIPTION 

Integrator with BCD Interface 

Microcomputer System with 
RS232 Interface 

Microcomputer ROM Expander 

Real Time Clock ROM Pack 

Binary /BCD I/O Interfaces 

- Interconnected 

- User Supplied Interface Box 

Minicomputer with 9 Track Tape Drive 
Interface Box with T.F.D Photo Count Display 

2 Channel - 12 Bit D/A Converter * 

12 Bit 16 Channel Data Acquisition System * 


QUANTITY MANUFACTURER & MODEL NUMBER 

4 Eppley Laboratory, Model 411-6140 

2 Tektronix , Inc . , Model 4051 

1 Tektronix, Inc., Model 4051E01 

1 Trans Era, Model 641-RTC 

5 Trans Era, Model 632 BCD 

with Options 1 and 2 

1 Digital Equipment, PDP 11/34 

1 Designed and Built by D. D. Venable 

& R. W. Blakey 

1 Trans Era, Model 620 DAC 

1 Trans Era, Model 652 ADC 


* These two devices were used in the design and testing stages. 
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APPENDIX V 

ABSTRACT - BETA KAPPA CHI PRESENTATION 

COMPUTER AUTOMATION OF A SOLAR RADIATION 
LABORATORY. R. W. Blakey * and D. D. Venable , 

Physics Department, Hampton Institute, 

Hampton, VA 23668 

A solar radiation measurement laboratory that includes two 
precision spectra pyranome ters , a precision infrared radiometer, a 
normal incidence pyrheliorneter and an all sky camera has been automated 
to allow direct computer acquisition of insolation data. Signals from 
the solar instruments are Integrated and displayed on five digit light- 
en'd, t ting-diode displays. The Integrated signals are transferred via a 
general purpose Interface card to a microcomputer. We have designed, 
constructed and implemented hardware and software configurations that 
permit data acquisition, storage, transfer, and display. System 
reliability tests have been performed and mean time between failure and 
system down time have been characterized. 
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APPENDIX VI 


ESTIMATION Of VIbSR UNCERTAINTY 


An uncertainty of 20%, in the calibration aid digitization of full- 
resolution visual imagery , was obtained using the following reasoning? 

Assume tliat the 5 x b image array used in estimating the spaceward 
reflectance of the Earth-atmosphere system is about 3/5 land (reflectance 
- 0.15) and 2/5 water (reflectance - 0.04) for the local solar measurement 
site geography using Table 7 in Muench (1981); the net random error, after 
multiplying the land aid ocean reflectivities by the weights just given, 
for this array is 9.8% for 1 mi x 1 mi resolution. Assuming that the 
inverse proportionality between the net random error and the nadir-point 
resolution holds between )? mi x )£ mi and 1 rni x 1. mi as it does between 
1 mi x 1 mi and 4 mi x 4 mi, this translates to a net random error of 19.6% 
for /z rni x /? mi. I he Pythagorean sum of this with the 5% systematic 
calibration errors for the variable model of Muench (1981) yields 20.2% 
■which should be rounded to 20%. 


